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Abstract 

We study in detail deep inelastic scattering in the 't Hooft model. We are able to an- 
alytically check current conservation and to obtain analytic expressions for the matrix 
elements with relative precision 0{1/Q'^) for 1 — x ^ This allows us to compute 

the electron-meson differential cross section and its moments with precision. For 

the former we find maximal violations of quark-hadron duality, as it is expected for a 
large analysis. For the latter we find violations of the operator product expansion at 
next-to-leading order in the 1/Q^ expansion. 

PACS numbers: 12.38.Aw, 12.39.St, ll.lO.Kk, ll.15.Pg 
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1 Introduction 

At its birth Quantum Chromodynamics (QCD) looked like a rather peculiar theory. It is 
constructed in terms of quarks and gluons, whereas all that one observes experimentally 
are hadrons, very specific combinations of those "elementary" degrees of freedom. Indeed, 
when the idea of quarks and gluons was first proposed [1], they were considered a mere 
fictitious tool to try to describe the hadron phenomenology. Nowadays, no one doubts 
their actual existence, as they leave their footprint in Deep Inelastic Scattering (DIS) 
experiments with hadrons, or in the ratio R = crg+e-^hadrons/ c"e+e-^^j+/x-; for example. Nor 
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does anyone doubt that QCD is the correct theory to explain their dynamics. However, 
thirty-six years after QCD was vindicated as the theory of strong interactions [2], we still 
lack a satisfactory analytic description of the hadrons in terms of the degrees of freedom 
and parameters that appear in its lagrangian. 

The difficulty resides in the fact that, leaving aside symmetry considerations (or how 
symmetries are realized), the only quantitative and analytic computational scheme to 
check the dynamics of QCD from first principles consists in weak-coupling computa- 
tions. In principle, those are limited to the computation of Green functions in the Deep 
Euclidean limit. The connection with experiment, however, requires the treatment of 
non-perturbative effects as well, and to relate those computations done in the Euclidean 
domain to the physical cut. 

Non-perturbative effects are taken into account through perturbative factorization 
techniques. The idea behind this approach is to try to separate the non-perturbative 
effects from the perturbative ones, dividing our calculations into two pieces: one which 
we can calculate perturbatively, and another which we leave unevaluated and determine 
through comparison with experiment or lattice calculations, for example. Essentially, 
all these factorization techniques are inspired on Wilson's Operator Product Expansion 
(OPE) [3]. The basis of the OPE is the application of the following relation in the deep 
Euclidean regioiil], 

I I e^^-dxT{A{x)B{Q)) ^ C„^^(g, ;,)0„(/x), (1) 

n 

where A, B are some local operators, 0„ are local operators with increasing dimensionality 
in n and the right quantum numbers to reproduce the left-hand side, C^^ are distribu- 
tions, and /i is the renormalization scale. The coefficients C^^ encode the physics beyond 
the scale /x, and the local operators encode the physics below this scale. In QCD the 
coefficients are calculated using perturbation theory, and the operators are assumed to 
hold all the non-perturbative physics. This is not completely accurate, however, as per- 
turbative effects can enter the matrix elements of the operators between the initial and 
final states, and non-perturbative effects make their way as well into the coefficients (for 
example, in the form of small-size instantons) . However, this is generally disregarded, 
and the OPE is used as a series of perturbative coefficients times some matrix elements 
to be determined experimentally or otherwise. The series is understood to be asymptotic: 
at some point in the expansion, non-perturbative effects are expected to cause it to break 
down [5]. Actually, the validity of the OPE is only established in perturbation theory [6]. 

^It should be mentioned that the primary definition of tlie OPE is without the time-ordering. The 
time-ordering introduces some ambiguities in the definition of the left-hand side of the equation (and 
consequently on the right-hand side). This is due to the fact that local terms in time are not fully 
determined (and we should also specify, in principle, in which frame we consider the time evolution) . As 
a matter of principle, one may try to fix them by asking the correlator to have the desired transformation 
properties under the symmetries of the system. In practice, we will consider the imaginary part of 
the correlator and obtain the complete result through dispersion relations. This guarantees the desired 
analytic properties for the correlator. 
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There is no mathematical proof that Eq. ([T]) can indeed reproduce well the (unknown) 
exact solution of QCD for processes which involve non-perturbative effects, even accepting 
its asymptotic nature. The validity of the OPE in these cases is just an assumption. 

The connection of the OPE with the physical cut can be performed through dispersion 
relations. This is a well defined procedure, and the sum rules obtained with it are as good 
as the OPE is. But the OPE, as stated above, is also used to directly compute quantities 
on the physical cut. If we knew the exact solution to QCD in the Deep Euclidean region 
(in a finite region) we could safely perform the analytic continuation from there to the 
physical cut, but as all we have at best are truncated expansions, this procedure can 
be a source of uncertainties, usually called quark-hadron duality violations (see [7] for a 
general discussion) . We stress that quark-hadron duality violations are usually disregarded 
without a good theoretical basis. Typically they are only discussed, sometimes, in analysis 
of the vacuum polarization, and even more scarcely in other processes like DIS or B decays 
like B Xg'j, see for instance pi^pHllTllTTjlT^ . Note that in these processes, perturbative 
factorization techniques, or the associated effective field theories like soft-collinear effective 
theory, simply neglect duality-violation effects completely. These effects can be easily seen 
in the large limit and quantified in the 't Hooft model (two dimensional QCD in the 
large limit [I3]). We do so here for the case of DIS (see [Tl] for the case of B decays). 

The practical version of the OPE (perturbative coefficients times non-perturbative op- 
erators I.5.J) is at the basis of computations at large Euclidean momentum of (the moments 
in) DIS and the vacuum polarization tensor, which so far have been thought to be among 
the more solid predictions of QCD, since they are not affected by quark-hadron duality 
problems. Therefore, the importance of setting the OPE and the factorization methods 
used in quantum field theories, especially in QCD, on solid theoretical ground can hardly 
be overemphasized. The OPE has been only partially checked in models, for instance in 
the 't Hooft model. This theory is superrenormalizable and asymptotically free, so it is 
a nice ground on which to test the OPEi. This was done at the lowest order in the OPE 
in Refs. [T51IT6] for the vacuum polarization and for DIS off a meson with nice agreement 
between the results of the model and the OPE expectations. In Ref. |17] the OPE was 
numerically checked in this model at next-to-leading order (NLO) in the 1/Q^ expansion, 
with logarithmic accuracy, for the vacuum polarization. In Ref. p!B] the main results for 
DIS at NLO were presented. In particular a violation of the OPE was found at NLO in 
the expansion. In this paper the details of that computation are presented. The 

paper is organized as follows. 

In section [2] we review the 't Hooft model. We will present the model, the semiclassical 
approximation to its solution [TOIITU] . and the transition matrix elements for a vector 
current (in two dimensions one can also obtain from them the matrix elements for the 
axial- vector current). 

In section |3] we study DIS in the 't Hooft model. We calculate the full, non-perturbative 
expression of the forward Compton scattering amplitude in terms of the 't Hooft wave 
functions and energies. As we mentioned, we observe maximal duality violations in the 

^In the 't Hooft model there are no marginal operators. Therefore, the coupling constant has dimen- 
sions and does not run; no renormalons should then arise. 
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physical cut when compared with the expression obtained from perturbative factorization. 
Analytic expressions for the matrix elements with l/Q^ precision for 1 — x ^ {3"^ IQ'^ are 
also given. We then compute the forward Compton tensor and expand it in the Deep 
Euclidean domain with l/Q^ precision. This result is compared to what we would obtain 
with the OPE. One would expect a perfect agreement at this order. However, surprisingly, 
we find that our expansion contains, besides the expected local matrix elements, some 
non-local ones at 0(1/Q^), which cannot be part of the OPE. These non-local matrix 
elements arise from the constructive interference between two (non- analytic) oscillating 
terms. 

In section m we present our conclusions. In appendix lAl we present corrections to some 
formulas of Ref. [H], where we studied duality violations in the context of semileptonic 
B decays in the 't Hooft model with I/toq precision. Nevertheless, the main conclusion 
of that paper remains unchanged. Namely, one observes no duality violations in the 
moments with I/tOq precision. 

2 QCDi^i in the large Nc limit 

The framework formed by QCD in two dimensions in the large Nc limit is usually called 
the 't Hooft model [13]. This model exhibits confinement: there are no free quarks, 
and the only states with finite mass are mesons (the mass of baryons grows with Nc), 
which are composed of exactly one quark and one antiquark, with an infinite ladder of 
gluons exchanged between them and an infinite series of "rainbow" radiative corrections 
to their propagators (in the large Nc limit only planar diagrams with no internal quark 
loops contribute, and in an appropriate gauge gluons don't interact with each other). 
In two space-time dimensions, the Dirac structure of the lagrangian becomes trivial and 
gluons can be integrated out, leaving us just with quark fields with no spinor structure. 
This allows us to solve the meson spectrum, which consists of an infinite tower of in- 
finitely narrow resonances, due to the large Nc limit, and features Regge behavior for 
large excitations. All of this makes the 't Hooft model an attractive framework where one 
can exactly "solve" QCD, and test computational techniques employed in the real world 
against exact results. 

In section I2TT] we will first present the appropriate coordinates and quantization frame 
to treat QCD in 1 + 1 dimensions, and in section [2^2] we will consider its large Nc limit, 
the 't Hooft model; in section 12.41 we will present the transition matrix elements for a 
vector current at leading order in (from which, in two dimensions, one can obtain 

the matrix elements for an axial- vector current). 

2.1 QCDi+i in the light front 

The QCD Lagrangian is given by 

^1+1 = -Ig%G''^^^ + i^YD, - m, + ze) , (2) 
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where = <9^ + ig^t^ and the index i labels the flavour. 

One usually works with Minkowskian coordinates, and quantizes the fields in the 
equal-time frame, where fields are defined at = constant. However, in some cases a 
different set of coordinates and a different quantization frame prove to be more useful. In 
the so-called light-cone coordinates, the Dirac structure of the lagrangian becomes trivial 
in two dimensions, and once everything is expressed in these coordinates, it is natural 
to choose a quantization frame in which fields are defined at a constant value of one of 
the light-cone coordinates, and not This is the light-cone quantization frame [20] . 
This quantization frame may be convenient when dealing with nearly massless particles. 
In four dimensions this line of research has been pursued by many groups, see pT] for a 
review. In two dimensions it can be seen that it is a natural framework on which to solve 
QCDi+i in the large Nc limit. 

2.1.1 Light-cone coordinates 

Let us define a basis in 1 + 1 dimensions with the two following light-like vectors (with 
the metric = g~~^ = 2 and zero elsewhere), 

n^ = (l,l), < = (1,-1). (3) 

Light-cone coordinates are defined like 

= ■ X = (a;° + x^) , x^ = ■ x = — x^) , (4) 



which implies that 



and 



X 



° ~ - (x^ + X ) , x^ = - [x^ — X ) , (5) 



d d d d 

d~ = 2—— = — + — = do + di ^ p- , d+ = 2j— = — - — = 9o - <9i ~ p+ , (6) 
ox^ ox^ ox^ ox ox^ ox^ 

P+x- P-x+ 

P ■x = \ , (7) 

2 2' ^ ' 

d^x = ]^dx+dx~d^~^x^ . (8) 

For the Dirac algebra it is useful to define the corresping light-cone matrices 

i+ = 7^, ?i- = 7" ■ (9) 

To have explicit expressions, it is useful to work with an explicit representation of the 
Dirac algebra. We will use the following Weyl-like representation for the Dirac algebra 
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so that the corresponding hght-cone matrices are given by 



7- = /-y = -2^(^° J) 7^ = 7° + y = 2z(^; j . (11) 
We can define as well the following projection operators (75 = 7°7^), 

2 2 4^+^" 4^ ^ vol 

If we write 

then the projection operators separate the two components of the field, 

A,^= (■*„-) , A-^=(;_). (14) 

2.1.2 The QCD lagrangian in the hght-cone frame 

Once the light-cone coordinates have been defined, the next step is choosing a certain 
gauge, the light-cone gauge. In light-cone coordinates, the gluonic field is represented by 
the components 

= n+ ■ A", A"'- = n-- A" . (15) 

The light-cone gauge consists in fixing A"'+(x) = 0; the reason for this choice will become 
evident in the next lines. In this gauge the QCD lagrangian in two dimensions can be 
written like 

£1+1 = i (5+A'''-)' + ^^l,+i/^-^.,++^l,_«5+^i,_-mi (^^|^+(-#i,_) + (-#,,_) , 

(16) 

where i is the flavour index. Now, quantizing in the light-cone frame consists in deflning 
the fields in this lagrangian at x'^ = constant. The coordinate plays therefore the 
role of time, the role of the energy being played by the conjugated variable P~. The 
other variables, P"*" (and P± in four dimensions) arc kinematical. For instance, the 
component of an hadron H behaves in a "free" -particle way, 

P^ = J2pt^ (17) 

i 

where the sum extends over all the partonic components of the bound state. This allows 
one to define the variable "xi" , which measures the fraction of P^ momentum carried by 
a given parton, 

x,^^. (18) 



In this quantization frame the field ip- is not dynamical (it doesn't evolve with "time" , 
and is therefore a constraint) and can be integrated out, 

TTl ■ 

In our gauge, gluons, represented by the component A~ , are non-dynamical and can be 
integrated out as well. After removing the constraints, the resulting Lagrangian can be 
written like 

i i 



where we have defined 



-1 , a; < , 

e{x) = { , x = 0, (21) 
1 , X > 0. 

Once we have the Lagrangian we can construct the Hamiltonian, 

~Y.T [ ^a;-rfrV^J,+t>.,+ (x-,x+)|x- -y-|V^J_+t>,-+(l/-,a;+) . (22) 

The representation of the quarks in terms of free fields in the light-cone quantization 
frame is (note that this assumes that > 0) 

M^) = / ^ {a{p)e-''' + b\p)e^'^) , (23) 

and the anticommuting relations are 

{a(p), at(g)} = {b{p), b^{q)} = 2{2ir)6{p+ - g+) , 

{a{p),b^{q)} = {b{p),a^{q)} = 0. (24) 



The free propagator in the hght-cone quantization frame looks hke 

/j2 1 U-\- 1~ j_ ''^i 7"*" j_ 



•^One should not forget that there is another constraint, the Gauss law, that restricts the Hilbert space 
of physical states to those which are singlet under gauge transformations. See for instance [22J, where 
one can also find a quantization in the path integral formulation. 
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Figure 1: The infinite series of "rainbow" radiative corrections to the quark propagator. 



where the / stands for free. The renormahzed propagator is given by the infinite sum 
of one-particle- irreducible diagrams. In the light-cone frame, and in the large Nc limit, 
these diagrams are limited to the rainbow-like diagrams shown in Figured! since no gluon 
lines can cross each other and there is no gluon self-interaction (the gluon lines in that 
diagram are not truly propagating in our gauge; strictly speaking, all gluon lines should 
begin and end at the same point in time). 

However, only the first diagram of this kind contributes. 




,27+ 1 



where we have defined (3"^ = g'^N(./{2TT). Adding a gluon line on top of this diagram 
produces a vanishing integral, which kills all the other "rainbow" diagrams. The infinite 
sum yields 



Recall that this expression is gauge-dependent. Should we have chosen to quantize in the 
equal-time frame, the expression of the renormalized propagator would be (again in the 
large A^^^ limit) 

p / xeg.time _ / « 1^ > 2 ^ k+ 2 ^ '"'^ /^oN 

' ~ J {2ny k+k- -mj + P^ + ie ' ^ ' 

The difference between the two propagators is (in momentum space) 

2 p 
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which illustrates the fact that the imaginary part of the propagator is independent of the 
quantization frame chosen. Note as well that this term is local in "time", proportional 
to 6{x~^). Such terms would jeopardize the expected covariance of the Green function. 
Let us note that if, for instance, we consider the OPE of such Green functions at leading 
order in we get something proportional to "^^q^/q^ in the equal-time quantization 

frame but ^~ lq~ in the light-cone quantization frame. This is not a problem by itself, 
since the propagator is not a physical quantity. 

2.2 The 't Hooft model 

By solving the eigenstate equation (taking into account the constraints, and using n to 
schematically label the quantum numbers of the bound state) 

p-\n) = p-\n), (30) 

one obtains the basis of states over which the Hilbert space of physical states can be 
spanned. Here we will focus on the meson sector of the Hilbert space and we will generi- 
cally label the state as n), where i labels the flavour of the valence quark, j labels the 
flavour of the valence antiquark and n labels the excitation of the bound state. 

The solution to Eq. fl30l) in the large limit gives us the spectrum in the 't Hooft 
model. In this limit the sectors with fixed number of quarks and antiquarks are conserved 
and consequently the number of mesons; in particular, the sector with only one meson is 
stable in the large iVc limit. Therefore, the bound state can be represented in the following 
way 

\^T.nT^ = 4f r ^^t'^" ffr) <aiP)bUPn - P)\vac) , (31) 
vNcJo V2(27r) \Pn/ 

where a is the color index, 0^ is a wave function representing the bound state, and the 
state is normalized as 

m\z'f; n)^ = (27r)2PW+<5„„5,,5,,v5(i-W+ - pW+) . (32) 

The superscript (0) stands for the large Nc limit, and Pji'^ the eigenvalue of n)^^'^ (we 
do not explicitely display the fiavour content of Pji'^ except in cases where it can produce 
confusion) . 

0^ can also be understood in terms of the gauge invariant "null-plane" matrix element 

^Hix) = J [ rfre^^"^"^(^ac|^,,+ (O,O)0(O,y-)^,,+ (r,O)K,j;n)W , (33) 

where $(x^,?/^) is a Wilson line, 

$(x-, y-) = P[e^''^y- '^--^+(--))] . (34) 
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P is a path-ordering operator. We have inserted the Wilson hne between the quark fields 
to make gauge invariance exphcit, although in the light-cone gauge (A"*" = 0) its expression 
is trivial. 

The fact that the number of particles is quasi-conserved could make it possible to 
formulate the theory along similar lines to those of pNRQCD (for a review see [23]). 
where the wave function (the 't Hooft wave function in our case) is promoted to the 
status of being the field representing the bound state. We will not pursue this line of 
research here, however. 

From the large Nc limit solution we can obtain the general solution to Eq. ( l30l) within 
a systematic expansion inl/Nc using standard "time" -independent quantum perturbation 
theory (we use as "time"). It has the following structure (the momentum of the bound 
state will not be displayed explicitly unless necessary) 

m,n' k Pm P^i 

+0 (±) . (35) 

where the second term in the expression is 1/ ^/Nl suppressed. Here we have used the fact 
that, at order l/y/TT^, P~ only connects neighboring sectors (n-mesons n± 1-mesons), 
becoming an almost diagonal infinite dimensional matrix (see also [21]). 



2.3 The 't Hooft equation 

By applying the operator P^ to its eigenstate \n) at leading order in l/Nc one obtains 
the 't Hooft equation 

^ + 7^ 0n^(^) - / dycl>^{y)F- ^ , (36) 

X 1-xJ Jo [y-xy 

where M„ is the bound state mass, x =_P^ /Pn i being the momentum of the quark i, 
and P stands for Cauchy's Principal Partcl. The renormalized mass is given by m^ji = rn^ — 
The principal value prescription serves to regulate the singularity of the integrand, 
which originates in the infrared divergence of the gluon propagator. 

This equation cannot be solved analytically in general, but much can be said about 
the wave function (t)n{x) and the spectrum. The 't Hooft wavefunctions are chosen to be 
real and normalized to unity 

'rfx0f(a^)C(^)='5n„^, (37) 



1 1"^ 



One can use the following representation of this distribution, 



{x - yf 2 



1 1 



{x — y + ie)^ [x — y — ie)^ 



11 



and they vanish at the boundaries with the asymptotic behavior 



0^^(x) = <a;^'(l + o(x)) , a;^0, (38) 
and similarly for x ^ 1 (changing i ^ j and x ^ 1 — x), where (3i is the solution of 

m^i + cot 7T(3i = , (39) 

which in the massless limit approximates to 

P^ = —-^ + o{m,^. (40) 
vr p 

The only case in which the analytic solution of the 't Hooft equation is known is the 
ground state with massless quarks. The solution in that case is (f)o{x) = 1, which means 
that 

limc^ = l. (41) 

In principle, there are several ways to obtain this result. One can work along the lines of 
Ref. [2S] to obtain an approximate Schrodinger-like equation, which can be approximately 
solved for the ground state. Another possibility to fix the value Cq is by matching the 
solution (f)o = 1 and the solution 0o = CqX^' in the region of overlap (the latter is valid for 
a; <^ 1, whereas it can be approximated to a constant, Cq, for values larger than e~^^^\ 
which is a very small quantity for small masses. Therefore, there is a region on which the 
constant solutions: "cq", and "1", overlap and should be equal by continuity). One can 
also use the value of limmi^o Jq 4>o{x) = 1; to fix Cq. 

The wave functions also obey the following very useful symmetry relations jl5]. 



C(a;) = (-l)"C(l-x) , (42) 



m, I ' dx^^ = (-l)"m, dx^^ . (43) 



2.3.1 Semiclassical solution of the 't Hooft equation 

For large n (high excitations) one can obtain approximate analytic expressions both for 
the meson wave functions and the spectrum through a semiclassical computation. In the 
interval 1/n < x < 1 — 1/n the WKB method gives the following solution for the wave 
function (valid up to 0{l/n)), 

0^^'(x) ~ V2 sm[{n + 1)'kx + 6ii{x)] , (44) 

and the spectrum reads [13] 



M'^ = n^P'n + {mU + ml^)\nn+^-^ + c(^]+c(^]+o(^] . (45) 



n 



12 



The phase shift 5n and the constant C{mjj^) in the spectrum were obtained in Ref. [19], 
and they are (the expression for C(m?^) was incorrectly written in Ref. 



siHx) 



1 

7T 



mi 



R 



[(1 — x)\n.n + Inx] 



X] 



C 



m. 



i,R 



2 



[x\nn + ln(l — x)] + x 



C 



m 



2 



TT 



TT- 



c 



P' 



m. 



i,R 



P' 



dy 



2y/smh.2y 



1 



ycothy + m^^/ P'^ y + n'^ 



+ 



^Ir 



In the hmit of small bare masses, this coefficient reads (C(— 1 

C 



P' ■ 
6.07242) 



^Ir 



P' 



C{-1) + v^TT^ + 0{m^ 



(46) 



(47) 



Eq. (jinD is obtained by studying the behavior of the wave function near the classical 
turning points, which is delicate for the 't Hooft model, and requires a precise quantum 
treatment of the boundary regions (0<a;< 1/ra, 0<1 — a;<l/r;,). This treatment is 
provided by the boundary- layer approximation. 



2.3.2 The boundary-layer approximation 

The boundary-layer approximation was first presented in Ref. [16] and later studied in 
Ref. [19]. As its name indicates, it is concerned with the behavior of the wave function 
on the boundaries, and it is valid only for large excitations. 
For X < 1/n the boundary- layer function is defined as 



(4J 



to 



for finite ^. For 1 — x < 1/n one may use the symmetry property shown in Eq. 
write 

0,(0= lim(-l)>^^Yl-e^l . (49) 



'M2 



0j(O approaches ^ = following the behavior of Eq. (l38l) . and matches the WKB solution 
(ilD as ^ ^ oo, 

0i(O^v^sin[(e/7r + 5XO] , if ^ ^ oo , (50) 

with 

i,K) = i|-!iiuKM-c(^V^^l. (51) 



r2 



The boundary-layer function fulfills the following equation: 



.(0 



"^IrIP 



(52) 
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It is possible to analytically solve the Mellin transform of this equation [19]. Define 

V;,(A) ^ / rfee'"Vi(0 • (53) 

In terms of ipi{X) the boundary- layer equation transforms into the following difference 
equation 

^i(A + 1) = (7rAcot(7rA) - Avrcot(7r/5i)) V^,(A) . (54) 

Any solution to this equation can have an arbitrary multiplicative periodic function 
P(A) = P(A + 1). A unique solution is determined by the following three conditions: 
i) analyticity of ip{XR + i\j) in the strip —[3 < A^j < 1, which can be seen by the definition 
of ip{\) and the behavior of 0i(O for small and large ^; ii) real analyticity, 'ipi^X*) = ^/'(A); 
and iii) the asymptotic behavior for Xj —>■ oo, which can be obtained from the WKB 
solution. With these conditions, the solution to Eq. fl5^ is 



p -"--V 1 + m: o/vrn^ tan7rA/( A + n) 



P " ' 1 + Kii/^/5')tan7rA/(A + n) 
where '?/'o(A) is the solution to Eq. (15^ in the massless case 



^o(A) = 7r"r(A)exp 



^ M + ^sin^vTM 
-27r / du 2 



sin(27ru) 



(56) 



and /5„ are the roots (0 ^ /3„ ^ 1) of 

7r/3^(/3„ + n)cot7r/3„ + m,^ = . (57) 
In particular, one obtains for A = 0, 1 

rrf«*«)=^f. (58) 

These two expressions are the leading order contributions to the integrals jj^ (ia;0^ (x)/x 
and dx(j)]{{x), respectively, for large n. Let us see that this is indeed the case. We 

consider the first integral; if we split it into two parts, cutting at some point /i such that 
1/n < fi < 1 — l/ra, we can write 

dx(j)'^{x)/x = / dx(f)]ii{x)/x + / dx(f)'^{x)/x 

Jo Jfi 

1 



dx(t)'^{x)/x + I dx(f)'^{x) 
dx(t^^{x)/x + {-lY I dxcpi^ix) 



'rf^M)+o(l/n), (59) 
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where in the third hne we have used the symmetry property given in Eq. ( H2l) . The integral 
dx(j)l{{x) can be obtained from dx(t)^^{x)/x through the 't Hooft equation. 
In Ref. [T7], by matching OPE and hadronic resuhs for two-point correlators, the 1/n 
corrections to the results of Eq. ( l58l) were found to be 



dx- 



\x 



X 



TT- 



[3 



rrii 



1 + 



+ 



-11 



O 



(60) 



and 



M't / dx(t)'^{x) 
'o 



1 + 



^Ir + ^Ir 



2n7r2/52 



+ 



nijm 



+ (-l)"7r/3mj- 



1 + 



^Ir + m 



2 



+ 



'-11 



+ 



(61) 



In Eq. fl60l) . the term ^^^(—1)" is the contribution from the boundary 1 — a; < 1/n, 
the term we discarded in the last line of Eq. (!59l) . The other correction is therefore 
purely the 1/n correction to the boundary-layer function. One should obtain and solve 
the boundary- layer equation at next-to-leading order in the 1/n expansion to obtain this 
correction analytically. Such a computation would require a dedicated study and goes 
beyond the aim of this work. 



2.4 Transition matrix elements 

In this paper we will need the transition matrix elements between mesons generated by 
the electromagnetic interaction, 

CQED (62) 

i 

This interaction does not change flavour. Therefore, we will only consider neutral cur- 
rents. We will consider the case of a charged meson made of a quark and antiquark with 
different flavour and the case of a neutral meson made of a quark and antiquark with 
the same flavour. The case of the charged meson is more interesting since it is stable 
under electromagnetic interactions. We obtain the matrix elements by using light-front 
Hamiltonian perturbation theory in the 1/Nc expansion, as we did in Ref. [14j for the 
case of the flavour- changing currents. For ease of reference we review the procedure here. 

We only aim to obtain the matrix elements at leading order in l/N^.. Nevertheless, 
this does not mean that we can just work with the leading-order solution to the bound 
states. As we will see, we will also need the l/\/W(. corrections to the bound state shown 
in Eq. (1351). 

The contribution to a matrix element can be split into two parts. We distinguish the 
contributions to the current according to whether they come from "diagonal" or "off- 
diagonal" terms, which we show in Fig. [21 The "diagonal" term directly connects the 
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Figure 2: Contributions to the hadronic matrix elements of the current (here shown for 
a general flavour- changing current ). The first figure corresponds to the "diagonal" con- 
tribution to the matrix element, Eq. [63^) . The second and third figures correspond to the 
"off- diagonal" term, Eq. ( [5^ . The ® represents the current, and the gluon exchange the 
effective four-fermion interaction in Eq. [2E) . 
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current to the leading 0{1/N^) term of the bound state. Considering a flavour- changing 
current for some generality, the "diagonal" term is 

{i'j]m\i!i'TiJi\ij;n)\diag. = ^^\i'j;m\4}i>T^i\ij;n)^^'^ . (63) 

This term is of 0{1/N^) and is produced from terms of the type ipi'Tipi ~ a|,ai + • ■ •, 
which in a way change the flavour of the quark from i to i' . Nevertheless, there is another 
possibility: ipi'Tipi ~ a|,6j (bi'tti), which can be understood as the creation (annihilation) 
of a new bound state. This possibility does not overlap with the leading-order term in the 
1/Nc expansion of the bound state, but it does overlap with the term. Whereas 

the matrix element connecting the one- meson sector with the two- meson sector is l/y/Wl 
suppressed, the overlap of the two-meson state with the current is ^/Nl enhanced. This 
is why this contribution has to be considered as well at leading order in l/Nc- We define 
the "off-diagonal" term as 

/dP^ 1 
2(27r)P+p(0)- P(°)- P(°)- 

X {vac\i>i>TilJi I a'; n') {n'; n' \ {i' j ] m\p- \ij ; n)^^K (64) 

This matrix element is the contribution of the term fej/Oj. It connects the 1 / y/Wl correction 
to the initial state (see Eq. (135|) ) with the leading final state. It is non-zero when 
P+ > P^. A bonus of working this way is that, once ^^\ii']n'\^^\i' j]m\P~\ii]n)^^^ 
has been computed, it can be used for any current. The other possible matrix element 
(involving the term a\,h\, therefore connecting the leading initial state with the l/v'iVc 
correction to the final state) contributes for P+ < P+. We will present here the matrix 
elements at leading order in l/Nc just for the case P^ > P^, which then read 

{i'i; m\ipiiTi)i\ij; n) = {i'j; m\^pi>Ti;i\ij; n)\diag. + {i'j; m\ipi'Tipi\ij; n)\off-diag. ■ (65) 

We can obtain the matrix elements for the case P+ < P+ simply by exchanging the labels 
n <->■ m in the expressions at the right-hand side of this equation. Nevertheless, we will 
not need them in the next sections. 

We present the matrix elements for a vector current, F = 7^^. From them one can also 
find the matrix elements for the axial- vector current, as this current can be expressed in 
two dimensions as a combination of the vector current and the tensor e'^": 

where = 1, or = 1. 
We define 

q = P^-Pn (67) 

X ^ -q^/p;l. 

The expressions below are therefore valid for x > 0. 
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(66) 



2.4.1 Neutral currents 



We will consider the case of a charged meson made of a quark and antiquark with different 
flavour and the case of a neutral meson made of a quark and antiquark with the same 
flavour. 



2.4.2 Charged meson: Non-equal mass case 

With our conventions, the full transition matrix elements for the "+" component of the 
current read 

t 



{ij;m\'ijja+ilji\ij;n) = 2{ij;m\'ijjl^'i/ji,+ \ij;n^ 



= 2P+(l-a:) 



dzcp^{z)cp^{x + {l-x)z) 



(68) 







x^P^ f f f dudvdz ,j}''^^_^''^'^l '^^}^,^ (O + (1 - x)z) - ci/^{xu)) 
Jo Jo Jo 



{x{i -u) + {i- x)zy 

{ij;m\i;j-f+i;j\ij; n) = -2P+(1 - x) dz(j)g{z)(pi^{z{l - x)) 



{<l>^^{z{l-x))-<i^^^{l-{l-u)x)) 



where 



Gu{u,v;q^) = 



n'=0 



- Ml 



(69) 



(70) 



For the "-" component we obtain 



2 

'P+ 



n 



(71) 







xl I dudz mx+il-x)z)-cp-ixu)) 



{ij;m\ijj-f i'j\ij]n) 



m 



/ /o 



{l-z){{l-z{l-x)) 



+ /3'(l-x) 



1 rl 



X 




dudz- 



^0 



C(^) [Q^Io dvGjjiu,v;q^) - 1 
{1 - z{l - x) - x{l - u)y 



mz{l-x))-c^^{l-{l-u)x)) 



(72) 
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where in Eqs. ( I7T!) and (172!) we have used that 



dv^^iv) = j dv(j)i;>{v) = , if n = odd , (73) 
and 

[dyf: M„^,^^^|^^ = q^ [ dyGuix, y; q') - 1 . (74) 

In order to obtain some of the above expressions we have also used the equations of motion 
in order to rewrite ipi- in terms of the physical component in the light-cone quantization 
frame, This is licit as far as the current is sandwiched between physical states. 
We note that the matrix elements are related by current conservation: 

q^{ij;m\tljh7~i^h\ij;n) = -q' {ij^mltphl^tphlij^n) , (75) 

which holds for arbitrary values of x and Q (and for any flavour). It is also useful 
sometimes to use —q^/q^ = Q"^ / {x'^{P^Y)). Eq. (1751) looks quite non-trivial if we take 
a look to the explicit expressions in Eqs. (1681) and (1711) . Nevertheless it can be shown 
to be an exact identity by a combined use of the identity Eq. ([71]) (this equality allows 
to rewrite the "off-diagonal" term in such a way that terms with a sum over infinity 
intermediate states drop out in the difference) and the 't Hooft equation, Eq. (136|) . 

Current conservation also implies that the vector current matrix element can be written 
in the following way 

{tr,m\^Pa''m:n) = ^P^ + PI^+ ^ "^ g^jP^;:(g^), (76) 

{tj;m\^P,Y^,\tj;n) = (^P„^ + P„^, + ^ " "^ ^J A^-^(g^) . (77) 

Obviously P^mil"^) ^nmi.'f') related by charge conjugation symmetry: 

j&f) = -{-ir^^p^:^{c]'). (78) 

This property can be easily visualized using the symmetry property of the 't Hooft func- 
tion given in Eq. (142 p . It allows us to easily write the antiparticle currents, Eqs. (1691) 
and (172 p . in terms of the particle currents, getting the correct —(—1)"-+™ factor. 

2.4.3 Chargeless meson: Equal mass case 

In the case where the particle and antiparticle component of the meson correspond to the 
same field, the expressions for the matrix elements can be simplified. By using 

{^^■m\i^^l'^^\n■.n) = (p.^ + + "2^"^ ?") ' . (79) 

we obtain 



19 



{n]m\'ilJi^^'ilji\ii]n) = 2(m; m|'?/'|, V^j.+ lu; n) 



[1 - 2i'+(l 



x] 



^0) 



-x'P^ f f f dudvdz ^ ^' [ct^^[x + (1 - - 

Jo Jo Jo ' ~ 



[1 



2 



^2 1 ' ^ j^{zmx+ii-x)z) 

" z(x + (1 - x)z) 



+ P\l-x) (81) 



-1 /.I 

X / / dudz- 
'o io 



(t>t{z) (q^ Jo dvGii{u, v; q^) - 1 



ixil - m) + (1 - x)z) 



\{x + {l-x)z)-4^:{xu)) 



The terms associated to "1" correspond to the particle current, and the ones associated 
to to the antiparticle one. Note that our results fulfill charge conjugation 

symmetry. 

The "+" component of the vector current was already computed in Ref. [16]. The 
computation of the rest of the matrix elements had to wait to Ref. p6], but we disagree 
with their results for the "-" component of the vector current. 



3 DIS in the 't Hooft model 

We consider here the differential cross section of the electron-meson scattering going to 
electron-l- anything: eM —>■ eX. The interaction we consider is therefore the one given in 
Eq. ([S2D- One should also have to include the leptons and photons, which we will not 
do explicitly. We will only consider the electromagnetic interaction perturbatively in e at 
the lowest non-trivial order. 

We are particularly interested in the situation when the momentum q transferred by 
the virtual photon is very large (DIS). Considering DIS with mesons rather than with 
baryons will allow us to use the results of the 't Hooft model. The meson has fiavour 
content M ~ qiqj. Unless explicitly stated, the formulas will hold true either if i = j or 
not. 

DIS in QCDi_|_i could be considered somewhat delicate, since Quantum Electrodynam- 
ics (QED) is confining in two dimensions. However, we will consider the electromagnetic 
interactions as pure current insertions. Working in the 't Hooft model we will be able 
to write down the full, non-perturbative expression of the scattering amplitude. As we 
mentioned, we observe maximal duality violations in the physical cut when compared 
with the expressions obtained from perturbative factorization. Analytical expressions for 
the matrix elements with l/Q^ precision for 1 — x ^ P'^/Q'^ are also given. We then 
compute the forward Compton tensor in the Deep Euclidean domain with l/Q^ precision 
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and compare the results to what we would obtain with the OPE. Surprisingly, we find 
that our expansion contains, besides the expected local matrix elements, some non-local 
ones at 0(1/(5^), which cannot be reproduced by the OPE. 

In sections 13.1.11 and 13 . 1 . 21 we will present the kinematics and definitions that we will 
use, and the expression of the scattering cross section; in section [3.1.3l we will give the ap- 
proximate form of the matrix elements with 1/Q^ precision; we will use these approximate 
matrix elements to give our result for the expansion of the forward Compton scattering 
amplitude in section 13.1.61 Next, in section 13.2.11 we will perform the perturbative cal- 
culation of the imaginary part of the amplitude; in section 13.2.21 we will obtain the OPE 
expansion from the previous result, and compare it with the exact expression. 

3.1 Hadronic computation 
3.1.1 Kinematics 




Figure 3: Deep-inelastic scattering off a light meson. The momentum of the photon q is 
incoming. 



The kinematics of DIS can be found in Fig. [31 They share some similarities with 
those of semileptonic B decays [H]. However, in B decays one has that = 0, which 
is an added constraint. Therefore, the expressions we have to deal with here are more 
complicated than those used in B decays, since now more kinematical freedom is allowed. 

The kinematical variables that we use are = —q^ > (this comes from the kine- 
matics of the scattering process) and (g = Pm — Pn) 

_ Q' _ g+ ...^ 

Xb = , x = -. (82) 

is explicitly Lorentz invariant, whereas x is the natural variable that appears in the 
solution of the 't Hooft equation. Momentum conservation in terms of x reads 

Pi = Q'^-^ + P^il - , (83) 

X 
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and in terms of xb reads 

Pl = Q'^-^ + Pl (84) 

Xb 



By combining both equations one obtains the following relation between x and xb 

X 



This equality is quite remarkable. It means that in two dimensions x is also Lorentz 
invariant. On the other hand, xb can be unambiguously obtained from Eq. ( !84l) : 

1 , , 

= p2_p2 • (86) 

With this expression we see that if the target is in the ground state (n = 0), is positive 
and runs from to 1. Having xb > means that we always have g° > 0, that is, the 
photon always transfers energy to the target. If the initial state is a resonance (stable 
under strong interactions in the 't Hooft model), we can have < P^, and the value 
of Xb is not restricted anymore. If P^ < Q^, xb will run from to some number greater 
than 1, but if P^ > Q^, xb can reach negative values (and therefore so does g°). We will 
not be concerned about this last situation since we will only consider the deep inelastic 
region, with Q'^ ^ P'^. 

The fact that xb can be larger than 1 may look surprising, compared with standard 
DIS in 4 dimensions. The reason has nothing to do with dimensions but with the fact 
that in four dimensions one usually considers the ground state as the initial state. In this 
situation the final state can only be the initial state or an excitation, which sets x^"^ = 1, 
whereas x'^°'^ > 1 happens for a final state with lower invariant mass than the initial 
state. 

If we now consider x, we see that Eq. ([H5D is quadratic in x, which means that for 
a given we have two possible solutions, x > and x < (that is, either g"*" < or 
g"*" > 0). This has a physical origin, the parity symmetry. The first solution corresponds 
to the frame where the initial lepton moves from right to left, and the second one to the 
frame where the lepton goes from left to right. We can give bounds to the values of 
and q~ in both cases through momentum conservation, and the requirement that P^ ^ 0. 
The conditions are symmetric, of course, 

< < 1 , ^ < ^ < oo left-moving incoming electron 

(87) 

% < < oo , < < 1 right-moving incoming electron . 

However, as we break parity symmetry in the light-cone quantization frame, the two 
cases will not be equivalent for us. With the transition matrix elements given in sec. 
12. 4[ it is easier to find approximate expressions for the cross section in the case of a left- 
moving incoming electron, and therefore this is the frame we choose. In our frame, then. 



22 



< xb < x^"^, where 

^r^ = — ^ > 1 , (88) 

and < X < 1, where the hmit x = 1 is only reached in the massless case when the initial 
and final state are in the ground state. 

Note that in Eq. (j83l) . as far as ^ P^, the factor P^{l — x) can always be neglected 
(in a first approximation) compared with Q'^^—^, independently of the value of x. This is 
not true in Eq. fl84p . where in the ^ 1 limit P^ can not be neglected. More generally, 
we can actually distinguish three different kinematical regimes: 



a) P^ = M^>Q' 

b) Pl = Ml^QA^cB 



c) Pl = Mlr^AlcD 



m 



> 



■mr^Q//3 , 1-xb~/9/Q, 
, l-XB</3yQ\ 

where m is the principal quantum number of the hadronic excitation. 

In this paper we choose to work in the target rest frame. The DIS limit corresponds 
to the limit where Q"^ ^ oo keeping xb fixed. In our frame this implies g"*" < 0, g~ ^ oo 
and the Bjorken xb goes like 

The Bjorken variable xb and x are in general different but they approach each other for 
large Q^. In general the momenta in this frame scale like 

P+ = P-~Aqc7D~/3 
q+ = -xP+ 

q- = Q^ixP^) 

P+ = P+(l-x) 

P- = P- + q- r^q- . (90) 

We have two possible expansion parameters, A = a/Aqcd/Q, which is always small, and 
A = \/l — X, which is small in the limit x ^ 1. We can see that in this limit the outgoing 
hadron behaves as a collinear particle, with a big momentum component P^ and a small 
invariant mass, although the scaling is not standard: 

P+ ~ A%cD , ~ QVAqcd . (91) 

There is another frame in which DIS is usually studied, the Breit frame, in which the 
photon carries no energy. The momentum components in this frame are 

g+ = -Q 

q- = Q 

= Q + ^^ 

Pn = I 
P+ = /+ 

Pm = Q + (92) 



23 



where are fixed by setting ~ Ql and ~ Ql'^, and 

In this frame the scahng is simpler: i+ ~ QA^ and ~ QA^, so that 

P+^Q\\P-^Q. (94) 

This is the natural frame to study DIS near x = 1 (there is no dependence on Aqcd 
neither in P+ nor in P^^, unlike in the target rest frame). The price is that here both the 
initial and final hadron are coUinear traveling in opposite directions. 

3.1.2 Scattering cross section 

The differential cross section is given by (///' represent the momentum of the incom- 
ing/outcoming lepton) 

(95) 

where S is the transition matrix operator, and \l) — Vl^ai{l)\vac). Expressing /'^ as 
I'ti — ;m _ can rewrite this as 

= MPn ■ P. - 2(t^fp^ 2(2.)^- I '-^l^l^- '}n^-mP.^^-P. 

(96) 

Being an observable quantity, a is gauge-independent, and we can choose the gauge for 
the electromagnetic field that we please. The usual choice is the gauge A~^j^ = 0, but 
there is also the option A'^j^ — 0. We can give compact expressions for a at LO in a in 
each of these gauges: 



a — 



(97) 



4(P„.p,-M>2)'/'y (27r)^ 
for the first gauge, and 



4(P„-p;-M>f)V2y (27r) 



(98) 

for the second one, where W^^ , l^^ , etc. are the components of the leptonic and hadronic 
tensors, defined as 

P'^ = le" J d^xe-''i-^{l\'iPi{x)Yi;i{x)'iPi{0)YMm , (99) 
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and 



^ h,h' 



m=0 

\2r2 



x(27r)V(P„ + g-P„) . (100) 

In the above sum over h and /i', there is only a contribution when both indices are equal 
to i, j. Since in the large Nc limit the spectrum is comprised of zero- width resonances, 
this tensor is a sum of deltas at the position of each resonance, a structure that cannot 
be reproduced by perturbation theory 

The leptonic tensor can be easily calculated, and at 0(q;°) we find 

Im[n = Im[ie2y"d2e-^«-^(Z|V;i(^)7>/(^)V5KO)7'^V'/(0)IO] (101) 

We can see that the imaginary part of the leptonic tensor obeys the identity 

^ Im[/++] = -4T^Im[r-] , (102) 



which implies 



^ W^^ = , (103) 



as required by charge conjugation symmetry. Current conservation in 1+1 dimensions 
also implies that the hadronic green function can be written in terms of one single scalar 
function (unlike in four dimensions, where we have two functions for a spin-zero particle) 

W'\g) = (p.- - (p.- - ^^^) W(g^ xs) . (104) 



Therefore, 



If one writes the hadronic form factor in term of the current matrix elements one 
obtains 

^ oo 

W{Q\xs) = -J2 \^^P^^i<l') + ejA^^i<l')\MiPn + q? - Pi). (106) 

m=0 
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where P^4' and A^^^ were defined in Eqs. (176|) and (1771) . 

The delta of momentum conservation imphes that W{Q'^,xb) 7^ only for xb > 
(because we set < Q^). We could define different hadronic tensors that would give the 
same result for the cross section, but would also have support for negative values of x^. 
Instead of defining W^" with J'^{x)J''{0) we could also define W^" using [J^(x), 7*^(0)], 
or W^" with {J^[x), J^(0)}. The scalar part of these alternative tensors would be 

W{Q^xb) = W{Q\xb)-W{Q\-xb) (107) 
WiQ^,XB) = WiQ^,XB) + WiQ^-XB) . (108) 

Note that W{Q^,xb) = -W{Q^,-xb) and W{Q^,xb) = W{Q^,-xb). Actually, the 
tensor W^" is often used to study DIS. 



3.1.3 Matrix elements in the DIS limit 

In section l2.4.1l we showed the matrix elements for fiavour- neutral currents. Our aim here 
is to obtain analytic expressions for these matrix elements with relative precision (9(1/(5^) 
in the situation when 3> (this means 1 — x > P/Q, and large m). In this situation 
we can use the boundary-layer function and its properties for the final state m. 

We first consider the "7"" current, for the more general case of a charged meson 
{i 7^ j). We will show only the particle matrix elements, as the antiparticle elements can 
be obtained from them using Eq. (ITHI) . The "diagonal" term is 

Pn Jo Z{X + (1 - X)z) 

We are only interested to compute these matrix elements for values of x = Xm satisfying 

= Q'^-—^ + Pl{l - Xm) » , (110) 

Xm 

since this is the requirement imposed by the delta of momentum conservation of W^^ . In 
this situation we can assume that ~ Q^^^^ and {z = ^ (3'^ /M'^) 

1 Xm ^ /? 1 Xm /-0 -I 

^ '^XfS" -^7yi^^- 

Xm ^^^ra Xm ^ 

Note however that this had not been true if we had kept m fixed but independent of x 
and we had performed the limit x ^ 0. 

We can then expand Eq. (11091) by considering that the wave function 4>^^{z) oscillates 
heavily except in a small region around the origin. In a boundary-layer-like fashion, we 
could say that the integrand is concentrated between and some finite ^ <C M^/ and 
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we can expand in z{l — <^ 1, 



PtJ^ \ J ~r m\ - J ' n \^ / 

^ \ X J Ci Z X 



2 



(1-x) / dzct>Z{z) + 



1 



(112) 



2 ^ ^2 1-x 



In the last two terms we have used that ~ Q 
The "off-diagonal" term is 

2 °° M^, 

{ij;n\ija-tPi\ij;n)\off-diag = -^P^i^ -^)Y1 ^2 ,1^2 (l^^) 



~'0 ^0 



{X{l-U) + (1 -X)2)2'"^" 



The main contributions to the integral over z must come from the end points of the 
function 05^(2;), as m is large and the wave function oscillates heavily anywhere else. We 
then focus on the regions 

m m 

with finite ^. The region z = 1 — is actually subleading. Then, considering that, 

in the sum over n', states with large n' are weighted more than those with a small n', we 
use the boundary layer for 0„/ (m) and cpn' (v) as well. We use Eq. ( 1741) to express the sum 
in terms of the Green function. Defining 

it can be shown that 



A/f2 /"^ m- r°° 



n'=0 

where 



(116) 



»i 

hlir])= lim ci^;Gii(r//3VQ^ t;, g^) . (117) 



^00 
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Combining all this information we can approximate the "off-diagonal" term by 
{i]]n\^a iJi\iJ]n)off-diag ^ -2 — T7^a;0„^(x)— / dt, dr] du- 



^ -2vr/5^^x0:^(x), (118) 

up to o (l/Q^) terms. The last equality is found assuming that the integral is dominated 
by the region ^,1/ 00, where we can approximate the behavior of the boundary- layer 

functions by \/2 sin((^/7r): 

. .00 .00 2sin(f/7r)sin(z//7r) . , , 

di dv dv I y' =^'- 

Jo Jo + + 

We have tried to confirm this result by numerically computing the integral 

f\J\.*&i^M^ , (120) 



Jo Jo Jo iz + x){z + y) 

which, for large n', tends to the integral in Eq. (1118^ . Unfortunately, the numerical 
computation of this integral is delicate (we use the method developed in [19] for our 
numerical calculations), as it requires a fine tuning between a small value of the integral 
and the large value of (for instance for n = 40 one has M^/f]"^ = 402.2). Moreover, 
for large values of n the integral becomes less precise. Nevertheless, the results we obtain 
appear to approximately converge (in an oscillating way) to the expected value, vr^, as 
can be seen in Fig. m 

Adding up the "diagonal" and the "off-diagonal" approximations, the total result for 
the "-" current reads (for x = Xm), 

— TTX ■ 

(zj; m\%lja~%lji\ii] n) = 27r/5- 



, , ^iR + mln ^ niinij ^ m,^ + (-l)'"m,m/ , 



+ imln + i-irm.rn,)^^':fixj 



1 



For the "-I-" current, Eq. ( ]68l) . we cannot proceed in a similar way since both the 
"diagonal" and "off-diagonal" terms are logarithmically divergent at 0{1/Q^). The "di- 
agonal" term in this case is 0(1/Q^), whereas the "off- diagonal" term is 0{1/Q^). It is 
possible to give an approximate expression for the diagonal term at the lowest non-trivial 
order 

2P+(l-x^)^ rf#j^(z)0;^(x^ + (l-a;^)z) =2a:^P+7r/?p0^^(x^) + o(^^^ , (122) 
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Figure 4: Numerical evaluation of Eq. U20\) ranging from n = 8 up to n = 40 . It oscillates 
around the expected value with increasing accuracy for increasing n until the numerical 
accuracy of the computation deteriorates. Calculations are done with rrii = rrij = (3. 



In order to reach the desired 0{1/Q^) accuracy, we rewrite the "+" matrix element such 
that the logarithmic behavior cancels (or in other words we rewrite it in terms of the "-" 
current using current conservation): 



{ij]m\ipa^i)i\ij]n) = x 



2(W 



{ij;m\ija iJi\ij]n) 



(123) 



which holds both for particle and antiparticle. 

Summarizing, we have obtained simplified analytic expressions for the "-" and "+" 
particle currents with relative accuracy 0(1/Q^) in the situation 1 — x ^ For both 

currents we can see that the off-diagonal term is a correction compared with the diagonal 
term. The antiparticle matrix element can be obtained from symmetry arguments. 



In order to have a complete control over the matrix element, we need an estimate for 
the region 1 — x ~ /J^/Q^ as well. When x — 1 we can approximate x ^ 1 — M^/Q"^, and 
express the matrix element for the "+" current as [16] 



{ij; m\i!i'y+ilJi\ij; n)\x~i-Ml/Q^,Q^^oc> 



Ml 



1+/3, r 



Ml JO 



dv 







(f + zf 



-hl{Mlv 



+ 



dz<j^{z){l - 

1 



Q 



2+2/3^ 



(124) 
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where 0^(a;) = c{(l — x)'^^ + o((l — x)^^). The matrix element in this hmit is suppressed 
by a relative factor l/Q"^^^ with respect to the leading term in Eq. (11231) (though it is 
enhanced with respect to the subleading term). Nevertheless, Eqs. (I123P and ( I124p cannot 
truly be compared since they refer to different regions in x. In any case, we will see that 
the associated contribution to the moments is subleading (as far as is not very large). 
This is due to the fact that it only contributes in a narrow portion of the total integral of 
the moment. 

In principle, Eq. (I123p and (11241) (at least the leading order expression) should merge 
in the intermediate region when 1 — x is small and yet ^ We can see that they 
do if we let increase in Eq. (11241) . 



3.1.4 The hadronic tensor 

The expression for the component of the hadronic tensor W reads 



W~ 



-y. 

rr. 



dP+ 



{ij;m\^eh^hiOh ^hiO)\ij;n) 



2(27r)P, 

{ij;m\^ehilJh{Oh'iph{0)\ij;n) 



(27r)V(g + P„-PJ 



SiMl-M^Sl-x)-Q'- 



,1 — 2; 



(125) 



Using Eq. (11211) we can give an approximate expression for W , valid for large and 
1 - X > /?VQ2, 



x) 



X 



1 H ^ TTT^ + 



2 



mi + (— l)™'mimj 



X 



ml^ + {-Vrm,m,)—A':^{\-x 



(126) 
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at 0(1/Q^), and consequently W reads 

/ /IT nnr \ ^ 

W 



x)-Q'- 



1 — a; 



X 



X < eiTUi 



a; 



-K^ + (-l)'"m,m,)— 0:^(1 -x) 



The expression in the equal mass case simphfies to 



(127) 



7 — ^;r^E^(<-^'(i 



-x)-Q' 



. 1 — X 



[1 _ (_!)"+-] <) ( 1 _ ^ + (1 + (_i)-) ^2 ^ 1 



Q 



2 



:i28i 



From W the component can be obtained immediately using Eq. fllOSp . 

At leading order in l/Q"^, approximating the mass of the bound state by ~ 
mvr^/?^, Eq. (11271) simplifies to (note that there is a discrepancy in the relative sign of the 
antiparticle contribution if we compare with the would-be analogous expression in Ref. 

I) 



W 



lead 



2x 



m ^ 



(129) 

At this stage one could try to approximate the sum over m by an integral by using the 
EulerMacLaurin formula at the leading order, i.e. replacing Ylm ~^ I This is neither 
mathematically justified nor can one quantify the error associated to this approximation, 
since we are dealing with Dirac deltas. For instance, it is not clear how to handle the 

interference term, which would go like ~ 0^(x)0:^*(x)(— 1)^^^^^ (though one could argue 
that it oscillates very quickly for Q"^ ^ 00 and away from the end points, the rate to 
which it vanishes can only be quantified when working with moments in the next section). 
Nevertheless, if we keep going and perform this naive averaging we obtain 



W 



LO 



^ ) " ( [e.m,0^^(x)] ' + {e,m,^t{\ - x)] ' ) , 



(130) 
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which agrees with the expression given by Einhorn [16]. Note however that its analytic 
structure is completely different from the one of Wiead- On the other hand this expression 
will be useful for us in the next sections. 



3.1.5 The forward Compton scattering amplitude 

In this section we will consider the amplitude for forward Compton scattering, which we 
will need for a comparison between our hadronic results and those from a calculation in 
perturbation theory. Let us define the tensor T^'^ in momentum space as 

T^'^iq) ^ t [ dV-^zr,n\T{f{x)f{0)}\zr,n) 



- [Pli - ^(Q^ xs) , (131) 

where in the second line we have used current conservation. Note that due to translational 
invariance 

T'^^(g) = T'^^'i-q) . (132) 

As the tensor structure is symmetric, this implies T{Q'^, xb) = T^Q"^, —xb)- By using the 
spectral decomposition of T'^" we obtain 

\mT{Q\ Xb) = 27iW{Q\ xb) = 2n(WiQ\ xb) + WiQ^, ~xb)) • (133) 

By using analyticity and the Cauchy theorem (and assuming that the contributions at 
infinity vanish fast enough), we can obtain the full functionality of T{Q'^,xb) from its 
imaginary part: let us define 

-P..-,/M.. = |^^. (134) 
then we can write {umin. = Q^/(2M„x^"^-)) 
TiQ^,XB) = 2 

. T^w-.f) . (135) 

where we have used that, since W{Q'^ ,xb) only when xb > 0, we only need the 
imaginary part on the positive branch of xb- This means that we could obtain the 
expression for T using only experimental data. 
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3.1.6 Moments at next-to-leading order 

Neither W{Q'^,xb) nor T{Q'^,xb) can be reproduced througli a calculation done in per- 
turbation theory for physical values of a;^ (0 < < x'b°'^'). W is a. sum of deltas, and its 
structure determines that of T, but perturbation theory yields a smooth function for W, 
as we will see shortly. A comparison between hadronic and perturbative results is only 
possible in the Deep Euclidean region, where both and xb are large and a perturbative 
calculation in QCD could be justified. 

T{Q'^,xb) admits an analytic expansion in 1/xb for xb > x^"^', 

T{Q',xb)=4: Y1 Mn{Q')^, (136) 

Af=0,2,4,... -B 

where Mn are the moments of W (we take this equality also as a definition for an arbitrary 



M^iQ^) = / dxBx'^-'WiQ^XB) 
Jo 



1 + x^-' — 
dx ^ ^ ^W{Q'\x) . (137) 



Note that only even powers of N appear in T, due to the symmetry property shown in 
Eq. (I132p . The above expressions can also be rewritten as an expansion around 1^ = 0, 



7V=0,2,4,., 



T{u,Q')=4 X ^^T(^HQ'), (138) 

where 



T^^'^Q") - r -^^^^^Q") = C-^-YMNiQ') . (139) 

Due to the structure of deltas in W, that of the moments will be a sum over m. We can 
find an approximate result for this sum through the Euler-Maclaurin expansion (taking 
52 = 1/6, 54 = -1/30,...), 

E /™ = / + 2 (/(O) + /K))+E 7^ if'-''-'\^l - f^''~'\0)) , (140) 

m=0 k=l ^ 

where /^"^ means the n-th derivative of the function. The limit m = corresponds to 
X = Xmax.i and the limit m*, when m* oo, corresponds to x = 0. At the limit x = 



(m* oo) we can use our expression for W given in Eq. (11271) . As 0^(x) ^^-^ c\x 
both f{oo) and /^"^(oo) go to zero. At the limit x = Xmax. = 0) we take Eq. fll24p to 
express the matrix element. With it we can see that both /(O) and /"-"^(O) are suppressed 
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by a relative factor 0{{[3'^ jQ'^Y^'^^^) with respect to the leading term, and so we discard 
them. The integral runs over all possible values of m, so in principle we should divide it 
into two regions, one in which the matrix element can be approximated by its boundary- 
layer expression (Eq. (11211) ). another in which the matrix element is given by Eq. (11241) . 
However, no matter whether we use one expression or another, the contribution from the 
end-point region is suppressed again by a factor O^^P"^ /Q'^Y^'^^^), due to the behavior of 
the wave function in that limit and the smallness of the region. So, at 0(1/Q^) we just 
change — >■ J dm in Eq. (I127p . and insert it inside the integral over x. At this order 
we can also use the asymptotic form of the spectrum for given in Eq. ( H5l) . With 
these approximations the expression for the moments is 



N 



dx 



X 



(C(^))' + e>?(0^^(l-x))^ 



m. 

Q 



mlj,d(Pli{x) 



dx 



-2e?m|0^(l - x) 



m 



»«(l-rc) 



dx 



1 1 



4-2p p ■ - ^ 



' 1 — 2r d 

2- 0tf (x)0^^(l -x)-x- i<l>^ix)<l^a - x)) 

1 — X dx 



(141) 

where the superscript NLO stands for "next-to-leading order" and means that this expres- 
sion is correct with relative 1/Q^ precision at finite A^. We have neglected the oscillating 
(— 1)™" terms, as they give a contribution suppressed by a relative factor of {P'^/Q'^Y^^^'^^^ . 
This is easy to see if one divides 



5^(-l)™/(m)= /M- E 



(142) 



m=0 



m,even 



m,odd 



and then applies the Euler-Maclaurin expansion to each separate sum. The leading con- 
tribution to the sums (the integrals) will cancel out, leaving only the sub leading ones. 
However, the product of two oscillating terms goes like (— l)^"* = 1, which gives rise to 
the interference term in the last line of Eq. (11411) . Being sub leading in l/Q^, this term was 
not considered in previous analysis P^P?] . We will see the importance of this interference 
term in the next section. 

It must be noted that Eq. (11411) is not valid for all values of A^. The factor in 
the definition of effectively selects the region of x that contributes the most to the 
integral. This is easily seen if we express x^ as 



X 



N 



^Nhi{l-{l-x)) 



-N{l-x)+0{{l~x)^) 



(143) 



As A^ — > oo, only the region 1 — x < 1/A^ will give a sizable contribution. As Eq. (11410 
assumes that the region 1 — x > P/Q dominates the integral, it is only vahd for N < Q/p. 
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To be more precise, for finite (though otherwise it could be large) the precision of our 
calculation is l/Q^, if scales with Q, the precision of our computation deteriorates, in 
particular for N Q/jS, the precision of our computation would be 1/Q, since there are 
(in principle) terms of 0{N (3'^ / Q^) ~ which we have not considered. 



NLO 



3.1.7 Determination of T 

Since we have approximate expressions for the moments from Eq. fll4ip . one may think 
that (at least an approximate expression for) T{Q'^,xb) could be recovered from them 
using Eq. ffT36D : 



T^^°(Q^XB)^4 Yl M^^^iQ')^. (144) 

AT—n 9 A *^-B 



N=0,2,4.. 

Nevertheless, this is not correct, or it rather should be quantified in which sense T^^'^ 
provides with a good approximation of T. Since, M^^^ is only valid for finite (but 
otherwise large), as compared with Eq. 01441) is only a good approximation of T 

for Q"^ and xb large. This means far away from the physical cut (T^^*^ is real in the 
real axis in this region). On the other hand T^^^ can be considered to be the generating 
functional for the moments with not very large A^. Moreover, it is useful to consider 
T^^*-* as a function in the xb complex plane by analytic continuation for the subsequent 
comparison with the computations using perturbative factorization. In this way T^^^ 
can be written in the following way 

T^"°(Q^x^) = 4 / ' rfyB-T^'''°(Q^yB) — — , (145) 

Jo VB 1 _ f^UR^ _ 

where H^'^^^ is given by Eq. (11271) . performing the substitution J2m ~^ I '^"^ {W^'^ is 
given in Eq. (11301) ). In this equality we have also fixed the behavior of the function in 
the physical cut (the imaginary part for real xb or the ie prescription) by demanding 
that it have the causality properties expected for a time-ordered propagator. Note that 
by approximating the sum over m by an integral we have lost the analytic structure of 
the imaginary part, since T^^*-* will have a continuous imaginary part, unlike that of the 
original T, which was a sum of deltas. However, as we have already mentioned, it is 
still interesting to consider the function T^^'^^ for in principle it should coincide with the 
result obtained from an OPE calculation, which we will perform in the next section. The 
function T^^'~' shares with the original T the trait that it is analytic everywhere on the 
complex plane except on the positive axis; therefore the resummation of the moments 
M^^° amounts to the computation of the dispersion relation. 

We can actually push the integration limits in Eq. (11451) to — oo and oo, respectively, 
in terms of the x variable. Since the 't Hooft functions cancel out of the interval (0, 1), 
all we are doing is extending the interval of integration over {xmax, 1)- By doing so, we 
are introducing an error of 0{1/Q'^y~^'^^\ which lies anyway beyond the accuracy of the 
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moments we are resumming. This way the comparison between T and T will be 
clearer. 

In order to give a compact, factorized expression for T^^*"* we defin^ the following 
functions (sums over the color indices of the fields are implicit): 



9.{y) = -±4>i{y) = {P^y 



dx 



r J 2(27r) 

X (^(n|V^J^_(x-)$(x-,0)^,,_(0)|n) + (n|V;l,_(0)V^,,_ 

2 



(146) 



m 



-im-y) = -{p: 



y 



dx 



^2 

P'^ x~ 



(0)|n) 



(0) 



-ly- 



2(27r) 

X ( (n|V^|_(0)<l>(0,x-)^,,_(x-)|n) + (n|^J__(0)^,v 



(147) 



n) 



(0) 



/oo 
n|^|__(a;")$(x", z")^/'^, _(2;")^j_+(0)<l>(0, z-)^/'i,+ (2;-)|zj; , 
-oo 

where as before the Wilson lines, trivial in the light-cone gauge, have been inserted to 
make gauge invariance explicit. The expressions for gi can be compared with those shown 
in Ref. [28]. We find agreement with them. 

The 13"^ /mf terms in Eqs. fll46p and (11471) have been inserted to cancel the "off- 
diagonal" contribution to the matrix elements (?T,|'?/'|_(a;~)$(a;~, 0)^/'i^_(0)|n) and 

(n|'i/'] _(0)$(0, x~)?/'j _(a;~) |n). Note however that this "off-diagonal" contribution does 
not correspond to the matrix element in Eq. (IMl) . The term hfaj in the current, from 
which the matrix element in Eq. (IMl) comes, gives a null contribution in this case. This 
is due to the insertion of an external momentum P^y. The "off-diagonal" contribution 
originates in this case from the term flj&j. The matrix element in Eq. (I148P does not 
have "off-diagonal" contributions: the "off- diagonal" contribution involves a two-meson 
intermediate state, as in Fig. [21 which is incompatible with the color structure of the 
matrix element. 

The functions g{y) encode non-perturbative information. We can write T^^'^ in terms 
of these functions times some other functions J which will hold the perturbative contri- 



'^Note that gi/j{y) = —^fi/jiy), where fi/j(jj) was defined in Ref. [THj. Note as well that the 
definitions of Ji are, accordingly, also slightly different here and in Ref. [18] . The difference between both 
definitions has to do on whether one chooses the "+" or "-" component for the distribution amplitudes. 
In particular, fi/j could be obtained from gi/j applying twice the equations of motion to their Fourier 
transform with the proper normalization. On the other hand the definition of gint. = fint. and Jint. are 
equal. 
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bution, 



T^^°{Q^,xb) = -2 (^-^^ j dy {e^iJi{x,y)gi{y) + epj{x,y)gj{y) + eiejJintXx,y)gintXy)} 

(149) 

where the functions J are defined as 



Jj{x,y) 
Jintix.y) 



X 



x^\l 



1714171^ 

2- ^ 



2 

2^ 

2 

2^ 



M2 „ 



+ X 



2'-^y' 



+ X 



Q"^ dx 



y 



2x\l-2y)-x^{l-y) 



Q'^ dx 
d_ 
dx 



— x^ + ie 

,3 



y 



y2 _ _^ 



— x^ + ie 



(150) 



This is the factorized form we expect from an OPE. However, as we will see shortly, the 
interference contribution (the term involving Qint.) lies beyond the domain of the OPE. 
The functions / are real, so the expression for ImT^'"^ has also a factorized form. 



ImT^L°(Q^x^) 



4 



dy {e^ilm[Ji{x,y)]gi{y) + epm[Jj{x,y)]gj{y) 
+eiejlm[Jint.ix,y)]gint.iy)} ■ 



Recall that we obtained the moments M^'"^, and therefore T^^'~', through the use of 
the dispersion relation given in Eq. fll35p . relating the discontinuity of T{Q'^, v) on the 
positive axis to its structure anywhere else on the complex plane. Therefore, the functions 
J are actually defined as (changing variables from J{x,y) to J{x, z/)) 



J{x, v) 



,„ ^ImJix, v') 



_ ^2 



:i52i 



le 



There is non-trivial information in this equation: the result from a direct calculation 
of J(x, z/) in the Euclidean and the result from the dispersion relation might differ in a 
polynomial (see e.g. [29]). The reason we use this prescription is that we are interested 
in the comparison between the hadronic and OPE results, and we will use the same 
prescription in the perturbative computation. In this way we aim to eliminate spurious 
differences between both computations due to dispersion-relation issues. 



Expression of M^^° in terms of matrix elements 

The expression for the moments given in Eq. fll4ip can be rewritten in terms of matrix 
elements, expectation values of some operators. We expect in this way to rewrite the 
moments in terms of an OPE expansion. However, not all the matrix elements that 
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will appear can correspond to an OPE expansion: the interference term can only be 
represented through non-local matrix elements, 

For simplicity's sake we will keep the factors of M^/Q^ explicit. They are actually 
not relevant for our purpose, which is spotting differences between full and OPE results. 
Although there is some dynamical, non-perturbative information encoded in M^, the 
presence of these factors in Eq. (11411) has a kinematical origin: they come from our 
definition of W (Eq. (11041) ) and the change of variables from xb^ox (Eq. fl85l) ) . Therefore, 
they will also be present in our later OPE calculation in exactly the same way. 

Starting either from Eq. (I14ip or from Eq. (11491) . using the definitions given in Eqs. 
(I146l) - (ll48p . and integrating by parts we can write M^^'^ in terms of matrix elements 
with derivatives inserted between the fields. The expression at 0(1/(5^) is 



QM '-^^'-^ ^ (P+) 



N+2 



2 



X 1 1 - 2^ - (iv + 2)^ + ^ W 1 ^-(0)1^^-^ 



e^(^j;n|V^[_(0)- 



-iD 



X I 1 - 2^ - (iV + 2)!fi + N^tlB^ I ^^.^_(0)|.,; n) 



-iD+\ -iD 
(iV + 4) 1-— — -2 



P+ ; p+ 

n I n 



x^0,z-)tP,,_{z-)^l40)^0,z-)ij,,4z~)\tj;n)^ 

(153) 

where ip^{x)(D+) = {D+'ip{x))\ and *D = |(z5 - D). 

The interference term is represented by the matrix element in the last two lines. The 
origin of the interference term is the constructive interference between two oscillating 
terms, (—1)™ ■ (—1)™ = (—1)^™ = 1. These oscillating terms give, after summing over m, 
an 0{j3'^ IQ'^Y^^^^^^ contribution to the moments, but their interference is enhanced. A 
non-analytic dependence on l/Q"^ like that of the oscillating terms seems out of the reach 
of an OPE. Note that for large m, one may think of 

(-1)™ e'^^ (154) 

which has a non-analytic expansion in 1/ j3. So, although the interference term is formally 
a simple NLO term in an l/Q"^ expansion, it is built out of terms which seem to be beyond 
an OPE expansion, whose non-OPE nature would survive in the form of a non-local matrix 
element. Either way, a complete understanding of this non-local 4-field correlator is still 
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lacking, but what is certain is that because of its non-local nature it is beyond an OPE 
expansion. 

The problem we are encountering might be enhanced by the large Nc limit: the (— I)*" 
terms arise from the null width of the resonances. Incorporating finite widths to the 
resonances (going to higher orders in the 1/Nc expansion), oscillations would be milder, 
perhaps moving this interference down to some higher order of l/Q^- In any case, our 
result seems to indicate a breakdown of the OPE for DIS in the 't Hooft model. 



3.2 Perturbative factorization 

In this section we will compute the amplitude for the forward Compton scattering fol- 
lowing the recipe of perturbative factorization. We will first compute the imaginary part 
of T , from which we will find the imaginary part of T*^^^, and then the full T^^^ 
through a dispersion relation. Perturbation theory only makes sense in the Deep Eu- 
clidean domain, but if T*^^^ is an analytic function, we can obtain its behavior in that 
region from its discontinuity in the positive axis, irrespective of whether this function 
properly describes T in the physical cut or not. 

We will compare our result for T^^^ with Eq. (11491) . From T'^^^ we will derive the 
coefficients M^^^ , which we will check against Eq. (I153p . 

3.2.1 Calculation of T'^^^ 
T is defined as 

h,h' ^ 

In perturbation theory, at leading order in P^, T reads 



OPE 



Ai I (fxe''^"'lef [/^+(x)(ij;n|^J__(x)$(x,0)V^i,_(0)|ij;n) 

+/^+(-x)(ij;n|V^|^_(0)$(0,x)^i,_(x)|ij;n) 

P+(-x)(zj; n|^] _(0)$(0,a;)Vi,-(a;)|u; n) 
+P^{x){i]] n|^j__(x)$(x, 0)^^- _(0)|ij; n) 

where P^{x) is the "+" component of the free quark propagator. 



P^ix) ^ -TrYi-pl{x)] 



(2^ 



2 



fc2 



(156) 



(157) 



We have written the Wilson lines to restore gauge invariance. These Wilson lines include 
both components of the gluon field, and A~ , in an obvious generalization of their 
definition in Eq. ([31 
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The picture for the process represented by Eq. (11561) is shown in Figure [5l for the 
particle case (the antiparticle case just involves switching around all the arrows in the 
quark lines). 




Figure 5: The direct and crossed diagrams representing, in momentum space, the two first 
terms in Eq. I^156\} . The momentum p is the momentum of the quark inside the meson n. 



In order to deal with the matrix elements we take advantage of the kinematics we have 
chosen: in our frame, q~ —>■ oo. Whether we consider the direct or the crossed diagram, 
this implies that the quark propagating between the two vertices has a very large p~ 
component, that is to say, its propagation takes place in a very short "time" (x"*" takes 
the role of time in our quantization frame). We can therefore expand the matrix elements 
in powers of x"*", 

{ij;n\ijl_{x)<!>{x,0)iJi,-{0)\ij;n) ~ (ij; n|^|_(x~, 0)$(x-, 0)V^i,_(0)|ij; n) (158) 

+ n\tljl_{x', 0)D-(!?{x-, 0)?/'»,-(0)|ij; n)x+ . 

Effectively, by letting the coordinate x~ untouched and expanding in x~^, what we are 
doing is resuming all powers of p~^ in an OPE expansion of the amplitude (recall that 
af^ ~ p+). Since we stay at 0{1/Q'^) we don't have to worry about the term in the 
second line in Eq. (I158P : due to the x~^ multiplying at the right its contribution will be 
suppressed by a relative factor 1/Q^. 

Using the definitions of Eqs. (I146p and (I147p we can reexpress the leading matrix 
elements as 

(zj;n|V^|,_(x-,0)$(x-,0)^,_(0)Kj;n) = ^ [ rfye^^%^^,(y)-^W(n|^l,_(0)V^,40)|n)(°) 

(159) 

(^j;n|4_(0)$(0,x-)^,,_(0,a;-)Kj;n) = l t/ye^^^^,(2/)-^(°)(n|^]__(0)^,- _(0)|n)W . 

(160) 

In writing ImTjpg, one could neglect the (3"^ /mf terms. The reason is that its Fourier 
transform is proportional to 5{y) and, after we sum the direct and crossed contributions, 
it will get multiplied by powers of giving a vanishing contribution to ImTjp^. Thus, 
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effectively, we are representing the matrix elements through some parton distribution 
functions Qijiy). We write the imaginary part of Tjp^ as 



ImTopg = ~(pTp j dygi{y)lm 



4e2 



oo 



+ ^|2 / 

OO 



dygjiy) im 



dygi{y) Im 



oo 



8e| '■^ 



— oo 




Next we consider corrections in Actually, we will write these corrections together 
with the leading order result we have just obtained in a combined single expression. We 
want to compute the following matrix element at 0{j3'^ /Q^): 

rf2xe'^-^(p|r{V^(x)7-V^(x)7/'(0)7-^/^(0)}|p) , (162) 

where \p) = a^'(j9)|0) (for the moment we do not specify the value of p^). We have to 
consider the diagrams in Figure [6l where the "blobs" represent the renormalized prop- 
agators to all orders in We only show the direct contribution, the crossed-diagram 
contributions can be obtained from it with the change x —x. A similar computation 
should also be carried out for the antiparticle contribution. Other possible corrections are 
given by the diagrams in Figure [7] (plus the symmetric ones), but these are suppressed by 
a relative order of 1/Q^. 
The first diagram gives 

nil mf 1 




(p+ + q+) {p- + q-) (p+)2 p+^q+_ + ,g 

^ ^fl + ^'j-^ , (163) 

where in the last line we have applied the momentum conservation delta to the factor 
"^i /((P"*" + Q'''')(P^ + 9 )) (to order /3^), since we are only interested in the imaginary part 
of this diagram. 
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Figure 6: Diagrams contributing to the perturbative compuation at 0{P /Q ) 



The vertex correction reads 




f 



(164) 



where we have defined 



' dy 1 

.y' (i_^)(i_^_^)_:|^^2 
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(165) 



X 



1 + 4^ 
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1-z 



Q 



2 • 



(166) 
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We then find 





{p+ + q+){p- + q-) {p+y p+ 



(167) 



^2 



2V m? QV (P+)%+ + „+ - + i£ ' 

where in the last line we have applied the momentum conservation delta coming from the 
intermediate propagator and expanded at order (this is the relevant result if we only 
want the imaginary part). Within this approximation we have m^/ {{p'^ +q~^){p~ +q~)) — 1 
and 

(1 + qVp'-)PVQ'f{-qVp^) - -J^ + + ^1^1 - fl + 2"^' 



(168) 

where in the last line we have only kept terms of order and approximated p^ ~ m^. 
It is interesting to discuss where this contribution comes from in the original integral in 
Eq. fll65p . Due to the delta of conservation 2; ~ 1 and one can rewrite z ~ 1 — 5. Then 
the integral has contributions from ?/ ~ 1 and ?/ ~ 5 -C 1. With the precision of our 
computation only the region y ^ S contributes to the integral (both regions would start 
to mix at order 1/Q^). 
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We can finally write the perturbative result for IniT . It reads 



ImT, 
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(169) 
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or (strictly at 0{f3yQ^)) 



dygjiy) Im 



2 \ 2 



+ ie 



ImT, 



8e 



OPE 



dygi{y) Im 



?/ 1-2 



8g2 ^oo 



{P^ 



+ ^2 



2^ 2A 



2^ 2A 



?/2 — ( 1 + 



?/2 — ^1 + 



,2x2 



(170) 



From ImTQp£,(Q2, x) we can obtain W oPEiQ"^ ^x), 



Wope{Q\x) 



2x 



-ImT^ 
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ehi {x {1+ ''^ 
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2 I / 1 , ' j,R 
^,9i (x(l + -^ 



Now, using dispersion relations (the analogous to Eq. f ll35p but with ±00 integration 
limits, since we do not need to know for which values of x the integrand gives a non-zero 
contribution, this is built-in in the result), expanding and integrating by parts, our result 
for r«^^ at 0(l/g2) is 



TO^^(Q^XB) 



dy {epi{x,y)gi{y) + epj{x,y)gj{y)} . (172) 



Eq. fll72p can be compared with our factorized expression for T^'"^ given in Eq. 
(11491) . We recognize the particle and antiparticle contributions, but the interference term 
is missing. 

The corrections involving the exchange of a gluon between particle and antiparticle 
represented in Figure [7] (plus the symmetric ones) are suppressed by a factor of 1/Q^. 
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Figure 8: Additional (non-perturbative) diagrams involving both the particle and the an- 
tiparticle. 

Perturbation theory does indeed contemplate interference terms, but they cannot account 
for gint.- 

Besides these perturbative corrections there are also the diagrams shown in Figure [H 
which are suppressed by a relative factor of and {P^/Q^Y'^^^ , respectively. 

Strictly speaking they do not belong to the dominion of perturbation theory, but they 
illustrate the fact that the interference term cannot be produced by any diagrammatic 
calculation at finite order in reinforcing the idea that our interference term, although 
formally an term, is non-OPE in nature. 

From our OPE calculation we see that we can understand the perturbative functions 
Ji and Jj as coming from the propagator of a collinear quark times kinematical and loop 
corrections (after an expansion). The connection between the functions J and the quark 
propagator that appears in Eq. (11721) is however not direct in general. Note, in particular, 
that the quark propagator may depend on the quantization frame and the gauge fixing 
used, which makes the explicit expression of the quark propagator different. We have 
tried to avoid as much as possible these ambiguities by using dispersion relations and 
demanding the T'^'^ to have the expected tensor structure. 



3.2.2 Moments 

From Eq. (11721) we can find the expression for the OPE coefficients M^^^. Actually, for 

this purpose it is more convenient to keep the factors of inside the functions fij as 
in Eq. (11711) . The expression for the moments obtained from a partonic approach is then 



Mr^(Q^ 



dx 
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(173) 
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which, using the definitions of gij given in Eqs. (I146P and (11471) and integrating by parts, 
can be reexpressed as 



gM '''' '^^'-^ ' (p+) 



7V+2 



(174) 



If we compare M^^'^ with M^^^ we confirm that the OPE does indeed get the terms 
involving local matrix elements right, but cannot grasp the non-local one. 



4 Conclusions 

We have thoroughly studied the 't Hooft model. We have obtained exact expressions for 
the current matrix elements in terms of the 't Hooft wave function. We have then studied 
the Deep Inelastic Scattering of a lepton off a meson in the 't Hooft model. We have 
calculated the full, non-perturbative expression of W'^'^ , and observed maximal duahty 
violations when compared with the expression obtained from perturbative factorization. 
Analytic expressions for the matrix elements with l/Q^ precision for ^ — x > (3/Q have 
also been given {1 — x > f3/Q means m 3> 1, where m is the principal quantum num- 
ber of the final hadronic state, so that we can use the boundary-layer function to find 
the approximate expressions for the matrix elements). This has allowed us to obtain 
expressions for the moments M/v at 0(1/ Q"^) for finite N (where we have also used the 
Euler-MacLaurin expansion). Here we have stumbled upon an unexpected result: the 
hadronic expression for the moments includes, besides the expected contributions from 
local matrix elements, a term at 0{1/Q'^) that can only be expressed through a non-local 
4-field correlator. This non-local matrix element represents the constructive interference 
of two oscillating terms, one from the particle and the other from the antiparticle. The 
oscillating terms go like (—1)™, and their contribution to the expansion is o (^^j) but 

their product doesn't oscillate and is of order 1/Q^. Upon resummation of the moments 
M^v we have found an approximate expression for T, T'^'"^, where T is the scalar part of 
the tensor T^'^. This result can also be obtained using dispersion relations, from W , 
which is obtained applying the Euler-MacLaurin expansion directly to W, the scalar part 
of W^'^ . T^^*-* is a good approximation to T for large and xb, and can be taken as a 
generating functional of the moments M^. 
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We have performed the same computation of W and T using perturbative factorization 
at one loop with precision. The perturbative calculation cannot see the structure of 
the bound states, and therefore VI^ope is a smooth function of x, unlike the hadronic result. 
A direct comparison between hadronic and perturbative results can only be performed 
in the Deep Euclidean region through the moments Mj^. From Wqpe we have obtained, 
through a dispersion relation, the coefficients of the OPE of T, the moments M^^^, 
with precision. We have checked that the OPE does get the contribution to the 

exact expansion from local matrix elements right, but misses the non-local one. This 
can also be seen at the level of dispersion relations, since W ^ W . Therefore, 
we conclude that this expansion breaks down at NLO for DIS in the 't Hooft model. 
The reason for this seems to be the non-analytic nature of the oscillating terms: their 
(— 1)™ behavior produces a contribution to the moments of 0{(3'^ / Q'^Y^^^^^^ , which seems 
out of the reach of an OPE. This non-OPE structure would survive in the (enhanced) 
interference term in the form of a non-local matrix element. We have considered diagrams 
representing interference between quark and antiquark, including some that are strictly 
out of perturbation theory's reach, and we have seen that neither of them could account 
for the interference term of the exact expansion, reinforcing the idea of its non-OPE 
nature. The acuteness of the problem might be due to the large Nc limit (therefore, 
one may suspect that one might run into the same difficulties in four-dimensional large 
Nc QCD). Key to the appearance of the interference term is the fact that resonances 
have zero width. With finite widths the behavior of the oscillations would be milder 
than (—I)™', which perhaps would move their interference down by some extra powers 
of l/Q^. Irrespective of this last comment, and in view of the findings of this paper, it 
is evident that more work should be devoted to a more rigorous study of quark-hadron 
duality and OPE-violation effects in perturbative factorization schemes. Otherwise the 
errors associated to those analysis will always have a certain degree of uncertainty, which, 
at present, cannot be quantified. 

The possible existence of OPE-breaking effects in QCD has already been discussed in 
the past. As early as in Ref. [5U] numerical evidence for the existence of OPE-breaking 
effects in the gluon condensate was claimed. Nevertheless, it is still unclear whether those 
effects can be associated to ultraviolet renormalons and/or higher orders in perturbation 
theory (for a recent discussion see [31]). Over the years there has also been some discussion 
on the possible existence of a {A^)min. condensate. This object should actually correspond 
to a non-local gauge-invariant condensate, though its explicit form is unknown for QCD 
|32] . Finally, there are some models that may produce effects that break the OPE, see for 
instance [33]. Nevertheless, those OPE-breaking effects would affect the static potential 
and the vacuum polarization. Regarding this we would like to emphasize that we do not 
find any OPE-breaking effect in the static potential or the vacuum polarization in the 't 
Hooft model. The static potential can be computed exactly in the 't Hooft model within 
perturbation theory. Therefore, there is no room there for effects associated to a sort of 
{A^)min. condensate. With the present precision of our computation, we also do not see 
OPE-breaking effects in the vacuum polarization [17]. Note that both in the case of the 
vacuum polarization and DIS we are talking about the same operator: the time-ordered 
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product of two currents. The difference comes from the physical states between which we 
sandwiched the operator: the vacuum in the first case and one particle state in the second. 
This may point to the fact that the OPE cannot be understood as an operator equality, 
as its validity may depend on the states between which the operators are sandwiched. 
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A Semileptonic B decays in the 't Hooft model 



In this appendix we present corrections to some of the formulas of Ref. pJ:J. There 
we studied duahty violations in the context of semileptonic B decays in the 't Hooft 
model with I/'hiq precision. The expression of the semileptonic differential decay rate was 
missing some terms, in particular the 1/n terms presented in Eqs. ( |60|) and (!6T!) . These 
corrections affect the computation of the moments but the main result remains unaltered, 
namely, one observes no duality violations in the moments with l/w^g precision. The 
reader is referred to Ref. [\A\ for definitions and a complete derivation, here we will only 
present the formulas and derivations that needed mending. 
In Ref. IHj we computed the decay rate 
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Approximate matrix elements 

We can expand the "diagonal" term of the matrix element of dV^^^ /dx for large m just as 
we did for DIS. If we define z = ^/?^/M^, the integrand will be concentrated on a small 
region of finite ^ near the origin. For the values of x such that z(l — x) ^ x, we can 
approximate 




Ml 

no 



dz- 



z)0g;(x + (l-a:)2:) 
zix + (1 — x)z) 



^ Ho 



X 



X) 



X 



dzcj)"^^ 



z - 



'He 



[X 



:i-x) 



x^ 



s,R _^ morris ^ 



xml + (— l)™'mcms 



M2 

m 



iQs 



M2 

m 



x) + O 



X 



M2 

m 




,Qs 
^Hn 



X) 



Ml 



(A.2) 

rnq. Note that in the last equality we could use 
Mhq{^ — x), since the physical matrix element is only defined for the values of x 



where we make the counting 



given by the delta of momentum conservation, 

Ml = Mlil-x^) , 



(A.3) 
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and/or use ~ mvr^/?^, since the above computation is meant for large values of m. 
Unlike in DIS, however, the condition z{l — x) ^ x does not always hold for physical 
values of x: in this case we have 

r M'^ r M'^ r ' ^ 

Thus, for Xjn ~ P'^/Mjj^, which corresponds to the largest excitations available for the 
decay, our approximation in Eq. ( ]A.2[) does not hold. In this region we should approximate 
the "diagonal" matrix element by 



^^^Hq Jo z 



(A.5) 



The contribution from this region is suppressed by the factor of , which goes to zero 
as niQ oo faster than l/mg, as Fig. [H] shows. Therefore, it can be neglected. 



0.004 



0.002 



, ^ 

20 30 40 50 60 ^ 

Figure 9: Plot of the evolution of c^^ with mq, compared with the curve given by l/niq, 

normalized to match the value of for mq = 10/3. The mass of the spectator antiquark 
is rus = p. 



The "off-diagonal" matrix element remains as it is shown in Ref. [Hj. Overall, we 
write the total= "diagonal" + "off-diagonal" matrix element in the following way for large 
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m and rriQ (up to a global sign, and setting x = Xm)'- 



IV 1 TT ^X 
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(A.6) 



We will need as well an approximate expression for the matrix element in the limit 
1 ^ ^ P'^/^Q (corresponding to the lowest resonances of the meson cs). Following the 
same procedure that led to Eq. fll24p in the case of DIS, we find 
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(A.7) 



This contribution is suppressed by a factor I/M"^^ with respect to the leading term in 
Eq. flA.6p . but it is enhanced with respect to the l/mg corrections, just like in DIS. 



For "intermediate" values of x (away from the boundaries x ~ jS'^/rriQ and 1 — x 
P'^/rriQ) the differential decay rate reads then, at ©(l/mg). 
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(A.8) 



Moments 



The differential decay rate is not a very well defined object in the large Nc-, since it becomes 
either infinity or zero. Its comparison with the expressions obtained from effective theories 
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that use perturbative factorization is not possible, as they yield a smooth function of x, 
so we turn to moments, which we define as: 
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(A.9) 



where Xmin. and Xmax. are given by 
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where Mm* is the mass of the maximum resonance of the meson cs allowed by momentum 
conservation, and Mq is the mass of its ground state. 

The calculation of the moments using Eq. flA.81) goes along the similar lines than our 
previous calculation of moments in DIS. As the differential decay rate is a series of deltas, 
the moments will be a sum over an index m that, as in DIS, we will rewrite using the 
Euler-Maclaurin expansion shown in Eq. (11401) . where now the limits m = m* and m = 



correspond to a; = Xmin. and x 



respectively. At the limit x 



Xr. 



(m 



m we 



can use the matrix element given in Eq. (1A.5|) . The contributions from both f{m*) and 
fi'^)[rn*) will be suppressed by a relative factor of 1/^q(c^q)^ < ^/^q with respect to 
the leading term; therefore, we neglect them. At x = Xmax. ijn = 0) we use Eq. (]A.7I) to 
represent the matrix element, and so we see that both /(O) and /'•"^(O) are suppressed by 
a relative factor of 0((/?^/mQ2)^"'"^^''). The integral goes from m = to m = m* (from 
X = Xmin. to X = Xmax.), SO principle we should divide it into three pieces: one for low 
values of X, in which we use Eq. (lA.Sp for the matrix element; another for intermediate 
values of x, in which we use Eq. (1A.6P : and another for high values of x, in which we use 
Eq. flA.7p . However, if we just insert Eq. flA.6P in the lower boundary, we will be making 
an error of o(l/mQ), and the contribution from the higher boundary is ©((/^^/mg)^"^^'^'') 
whether we use Eq. (IA.6|) or Eq. flA.7l) . 

Summing up, we take Eq. flA.Sp . make the change J2m ~^ I insert it into Eq. 

flA.gp . and what we obtain is (using Eq. fHSl) for the spectrum) 
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There has been a number of approximations here. In the first hne we neglected the 
oscillating (—1)™ terms, since their contribution is suppressed by a relative factor of 
o{(3'^ / w^q) . In the second line we have reshuffled the NLO correction in a way that is 

correct at the accuracy of the calculation. And in the last line we have used — 

^^Q,R '"■Q,R 

which is correct again with the accuracy of our calculation (this approximation is wrong 
as a; ^ 0, but there the asymptotic behavior of ensures that the error be of 

mm. I rn'^ 



-Ir 



0, but there the asymptotic behavior of cpfj^ 
o(l/mQ)), and we have also extended the lower limit of integration from x = Xmin. + 



too; 



1; in the lower limit the error 



2 N l+2l3s 
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to X = 0, and the upper limit from x = Xmax. + 

will be o(l/mQ), and in the upper limit it will be of 0((/?^/ 

The right-hand side of Eq. (lA.lip contains some implicit dependence on the heavy 
quark mass, since so far we have used the exact Hq-meson. If we perform an explicit 
expansion in l/mq, one obtains for the first moments, up to ©(l/mg). 
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(A.12) 

(A.13) 
(A. 14) 



where the static limit expectation values are defined in Ref. [14j. 

Eq. flA.lip is exactly the result we showed in Ref. [13] ■ Therefore, the conclusion 
we reached there for the moments still holds: they show no duality violations with I/toq 
precision. 
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